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ABSTRACT 

Magnetic Doppler imaging is currently the most powerful method of interpreting high- 
resolution spectropolarimetric observations of stars. This technique has provided the very 
first maps of stellar magnetic field topologies reconstructed from timeseries of full Stokes 
vector spectra, revealing the presence of small-scale magnetic fields on the surfaces of Ap 
stars. These studies were recently criticisied by Stift et al., who claimed that magnetic inver- 
sions are not robust and are seriously undermined by neglecting a feedback on the Stokes line 
profiles from the local atmospheric structure in the regions of enhanced metal abundance. We 
show that Stift et al. misinterpreted published magnetic Doppler imaging results and consis- 
tently neglected some of the most fundamental principles behind magnetic mapping. Using 
state-of-the-art opacity sampling model atmosphere and polarised radiative transfer codes, we 
demonstrate that the variation of atmospheric structure across the surface of a star with chemi- 
cal spots affects the local continuum intensity but is negligible for the normalised local Stokes 
profiles except for the rare situation of a very strong line in an extremely Fe-rich atmosphere. 
For the disk-integrated spectra of an Ap star with extreme abundance variations, we find that 
the assumption of a mean model atmosphere leads to moderate errors in Stokes I but is neg- 
ligible for the circular and linear polarisation spectra. Employing a new magnetic inversion 
code, which incorporates the horizontal variation of atmospheric structure induced by chem- 
ical spots, we reconstructed new maps of magnetic field and Fe abundance for the bright Ap 
star a 2 CVn. The resulting distribution of chemical spots changes insignificantly compared 
to the previous modelling based on a single model atmosphere, while the magnetic field ge- 
ometry does not change at all. This shows that the assertions by Stift et al. are exaggerated 
as a consequence of unreasonable assumptions and extrapolations, as well as methodologi- 
cal flaws and inconsistencies of their analysis. Our discussion proves that published magnetic 
inversions based on a mean stellar atmosphere are highly robust and reliable, and that the 
presence of small-scale magnetic field structures on the surfaces of Ap stars is indeed real. In- 
corporating horizontal variations of atmospheric structure in Doppler imaging can marginally 
improve reconstruction of abundance distributions for stars showing very large iron overabun- 
dances. But this costly technique is unnecessary for magnetic mapping with high-resolution 
polarisation spectra. 

Key words: stars: chemically peculiar - stars: atmospheres - stars: magnetic fields - stars: 
spots - stars: individual: a 2 CVn 



1 INTRODUCTION 

Doppler Imaging (DI) of stellar surfaces represents one of the 
spectacular successes of modern observational stellar astrophysics. 
Since the first applications of these tomographic modeling tech- 
niques to cool and hot stars in the 1970s and early 1980s 
(iKhokhlova & Rvabchikoval 1 19751 : iGoncharskii et all 1 1971 Il982l : 
IVogt & Penrod 19831) . they have revealed the complex surface 
structures of rotating stars in remarkable detail, providing impor- 



tant new insights into the processes occurring i n stellar outer layers 
and un iq ue constraints on t heoretical models dDo nati et al. 2003, 
2006 35 IPetit et all 120081: iMorin et alJ 120081 : iKochukhov etaO 
boollKochukhov & WadehoiCh . 

Conventional DI interprets rotational modulation of the in- 
tensity spectra of stars to reconstruct two dimensional maps of 
photospheric temperature or abundance spots (in the case of cool 
and hot stars, respectively). Basically, the stellar surface is divided 
into a grid of pixels, each of which is assigned independent lo- 
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cal physical properties (temperature or abundance). Synthetic disc- 
integrated spectra are computed, taking into account the local tem- 
perature/abundance of each pixel, as well as its rotational Doppler 
shift, and compared to a timeseries of observed spectral line pro- 
files. Iterative adjustment of the local physical conditions leads to a 
model (or "map") that is capable of reproducing the observed line 
profile shapes and modulation in detail. 

As the available observational data and computational re- 
sources have improved, DI models have become increasingly so- 
phisticated, with better resolution and realism (e.g. inclusion of 
various physical effects such as differential rotation). A qual- 
itative s tep forward was taken with the inclusion of m agnetic 
fields by IPiskunov & Khokhloval dl983h for Ap stars and bv lSemell 
dl989h for cool activ e stars. Zeeman or Magnetic Doppler Imag- 
ing (ZDI or MDI: iBrown et all Il99ll : IPiskunov & Kochukhov! 
120021 : Kochukhov & Pisk unovl 1200217 extends the basic princi- 
ples of DI to allow mapping of stellar magnetic fields by inter- 
preting timeseries of line profile Z eeman polarisation. Currently , 
the most sop histicated MDI co des (Piskunov & Kochukhov 2002; 
Kochukhov & Piskunov 2002) are capable of simultaneous and 
self-consistent mapping of surface abundances or temperature dis- 
tributions along with the vector magnetic field from timeseries of 
line profiles in two {IV) or four {IQUV) Sto kes parameters. 

Recently, IStift Leone & Cowlevl d2012h reported the results 
of experiments that seriously question the basic reliability of the 
results of DI and MDI, particular as applied to chemically pecu- 
liar Ap stars. As pointed out by Stift et al., DI techniques have 
revealed some surprising and exotic characteristics of the atmo- 
spheres of these objects. One of the most outstanding is the large 
chemical abundance contrasts in atmospheres of Ap stars, which 
according to DI analyses can be quite extreme (e.g. variations of 
several orders of magnitude across the stellar surface). In regions 
of high enrichment, some metals are sometimes inferred to be 
only a factor of ^ 3 less abundant compared to hydrogen (e.g., 
iKuschnig etall ll998). which is challenging to explain theoreti- 
cally. The distributions o f abundances can be extremely complex 
(iKochukhov e~a l. 2004b), while sometimes they are rather simple 
(e.gXuftinger et al. 20K)S). The distributions may show straight- 
forward relationships with the in ferred magnetic field geome- 
try (e.g. jRice. Wehlau & Holmgren! 1 19971) . but frequently they do 
not. Recently, MDI mapping using four Stokes param eter datasets 
(IKochukhov et a l. 2004a; Kochukhov & Wade 2010b has revealed 
the unexpected presence of intense, small-scale structures in the 
reconstructed magnetic field. The existence of these structures is 
unexplained, and adds a qualitatively new element to our growing 
picture of the magnetic stru cture of Ap stars. 

The essential claim of IStift et al.l j20ll hereafter SI 2) is that 
maps computed for Ap stars using DI and MDI methods are funda- 
mentally unreliable and that the variety of unexplained properties 
of maps described above are artefacts of neglecting the impact of 
large abundances and a bundance contrasts on t he local stellar atmo- 
spheric structure (e.g., lKhan& Shulvak l2007h . In particular, Stift 
et al. state that ignoring these effects results in artificial line profile 
variability, which translates into spurious abundance and magnetic 
structures in DI maps. Ultimately, they conclude that complex DI 
maps with high-contrast abundance distributions are incorrect, and 
that the presence of small-scale structure and high-contrast mag- 
netic spots, as inferred from MDI, cannot at present be reliably as- 
certained, and are likely an artefact of the inversion procedure. 

In this paper we confront the assumptions and methodology 
of SI 2, demonstrating that their conclusions are largely founded 
on fundamental errors and unrealistic assumptions, and conclud- 



ing that DI maps based on single, mean stellar atmospheres are 
highly robust and reliable. Our paper is organised as follows. First, 
in Section [2] we bring to the reader's attention several key aspects 
of DI and MDI methodology overlooked or misinterpreted by SI 2. 
Calculation of realistic opacity sampling model atmospheres, cor- 
responding local continuum intensities, and the Stokes IQUV pro- 
files is described in Section[3] Using the Fe abundance distribution 
and magnetic field topology reconstructed for a 2 CVn by KW10, 
we compute in Section |4~2l the full Stokes profiles of a strong Fe II 
line, exploring the consequences of the lateral variation of atmo- 
spheric structure on the intensity and polarisation profiles. A com- 
plete MDI inversion based on a grid of model atmospheres with 
different elemental abundances is presented for a 2 CVn in Sec- 
tion [43] We summarise and discuss results of our investigation in 
Section [5] 



2 GENERAL CONSIDERATIONS 

2.1 Robustness and stability of magnetic inversions 

As with any complex data modelling technique, capabilities 
and intrinsic limitations of MDI must be tested comprehen- 
sively. The physical foundations and numerical methods im- 
plemented in the In VERS 10 code were described in detail by 
IPiskunov & Kochukhov! (l200l hereafter PK02). This paper pre- 
sented a general discussion of the full Stokes vector inversion 
methodology and carried out a comprehensive theoretical and nu- 
merical assessment of its key algorithms, such as polarised ra- 
diative transfer and regularisation of the ill-posed stellar surface 
mapping inverse problems. The forward Stokes parameter profiles 
computed with InversIO were compared with the spectra pro- 
duced by two othe r independent magnetic spectrum synthesis codes 
dWade et al.ll200ll) . This study showed a very good agreement of 
the local and disk-integrated Stokes profiles produced by all three 
codes, provided they use a consistent definition of Stokes parame- 
ters and the same input mod el atmosphere and atomic data. 

A subsequent study by IKochukhov & Piskunovl d2002l here- 
after KP02) presented extensive numerical experiments designed to 
evaluate performance of Invers 10. These tests demonstrated that, 
given high-resolution IQUV observations, the magnetic inversion 
code is capable of correctly reconstructing abundance and mag- 
netic field vector surface distributions simultaneously and without 
any prior assumptions about the large-scale magnetic field geome- 
try. At the same time, it was found that an accurate mapping of a 
globally-organised magnetic field from circular polarisation alone 
requires an additional multipolar constraint on the possible mag- 
netic field structure. KP02 showed that the full Stokes vector-based 
simultaneous magnetic and abundance inversions are generally sta- 
ble and robust. It was found that the derived solutions are unique 
and the magnetic mapping code cannot arrive to a spurious solu- 
tion which would provide a good fit to observations. Another im- 
portant outcome of the MDI tests carried out by KP02 is that mag- 
netic inversions extract information about magnetic field primarily 
from the rotational modulation of Stokes profiles rather than from 
the Doppler shifts, unlike for the conventional scalar DI, making it 
possible to apply MDI to very slowly rotating stars. 

The comprehensive discussion of the MDI methodology im- 
plemented in InversIO and extensive numerical tests of this code 
presented by PK02 and KP02 are neglected by SI 2, who do not 
mention either of these two papers. This oversight makes much of 
the general criticism by S12 irrelevant because many of their con- 
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cerns about magnetic inversion methodology have been explicitly 
addressed by PK02 and KP02. 

We note that the contentions of S12 about what DI inver- 
sion can or cannot do lack a reference to the basic numerical tests 
demonstrating robustness of their own abundance inversion code. 
For instance, comparing Fig. 7 and 8 of their paper, one can see that 
the authors have failed to recover the true distribution of abundance 
spots even when they used a constant magnetic field and adopted 
the same mean atmospheric structure for calculations of the input 
spectra and for the inversion. A major discrepancy between the in- 
put and reconstructed abundance maps revealed in such a simple 
test indicates serious problems with the DI code used by SI 2. Ev- 
idently, results based on the application of this untested inversion 
code must be considered with caution. 

2.2 Regularisation in magnetic DI 

Typically, a grid of several hundred to several thousand surface el- 
ements is employed in DI. Because even a high-quality spectro- 
scopic data set may only represent a few hundred independent con- 
straints on the stellar surface structure, the conventional tempera- 
ture or abundance Doppler imaging problem is mathematical ly ill- 
posed dGoncharskii et alJll977l : IVoet. Penrod & Hatzeslll987l) . As 
a consequence, a large family of maps - each with the same large- 
scale structure, but differing in smaller-scale details - is able to 
reproduce the observations. In fact, the surface grid is often suf- 
ficiently fine as to be able to fit the observational noise. To solve 
these problems a penalty (or regularisation) function is included in 
the optimisation procedure, serving to limit the information con- 
tent (maximum entropy regularisation) or high-frequency structure 
(Tikhonov regularisation) of the map. It has been demonstrated that 
in conventional DI, for high resolution and signal-to-noise ratio 
data sets, t hese two approaches le ad to maps that are essentially 
equivalent ( Stra ssmeier et al.lll99lh . 

The situation is different for magnetic DI. IPiskunovl (2005) 
showed that the underlying inverse problem based on full Stokes 
observations is well-posed and thus has a unique solution, and is 
therefore not influenced by the selection of a specific type of reg- 
ularisation. In this case the role of regularisation is limited to sup- 
pressing numerical instabilities induced by a sparse wavelength and 
phase sampling typical of real observational data. In contrast, the 
MDI problem based on only Stokes I and V is intrinsically ill- 
posed, leading to different results depending on the applied regu- 
larisation algorithm. For instance, iBrown et all dl99lh have failed 
to recover a dipolar magnetic field distribution from Stokes IV data 
with their maximum entropy ZDI code but succeeded in recon- 
structing topologically more complex field geometries consisting 
of isolated spots. 

This counterintuitive situation when mapping simpler global 
field topologies from Stokes IV spectra is seemingly more diffi- 
cult than recoverin g very complex fields of late-type stars (e.g., 
Donati et al has been addressed by KP02. They showed 

that the maximum entropy constraint is inappropriate for global 
magnetic topologies since they do not correspond to a map with 
a lower information content compared to a distribution of iso- 
lated spots on a homogeneous background. KP02 achieved some- 
what better results with IV imaging of dipolar fields using the 
Tikhonov regularisation and demonstrated that a fully satisfac- 
tory reconstruction of the global field topologies is possible with 
multipolar regularisation, which penalises a deviation of the mag- 
netic geometry from a low-order multipolar solution. A simi- 
lar magnetic inversion technique, employing a general spheri- 



cal ha rmonic expansion, was independently developed by iDonatil 
(1200 lh and successfully used for the Sto kes IV ZDI of global 
fields of Ap stars f Luftinger et al. 2010 a]), fully convective low- 
mass stars dDonati et al.ll2006al : lMorin et alj2008D . s olar analogues 
(Petit et al. 2008), and late-type Ap-star descendants dAuriere et al.l 
12011 ). These developments were entirely neglected b y SI 2, who 
subscribed to the obsolete view of iBrownet al]([l991) that global 
magnetic fields can be only mapped using full Stokes data. 

Discussing the IQUV magnetic mapping of a 2 CVn by 
KW10, S12 insisted, citing Figs. 6 and 8 of KW10, that these in- 
versions are unstable because "results depend too strongly on the 
regularisation parameter" '. But these figures show exactly the op- 
posite, demonstrating the impact of changing the regularisation pa- 
rameter (which controls the relative contributions of the penalty 
function and chi- square of the fit) by a staggering factor of 10. This 
large increase of regularisation was required to wipe out the small- 
scale field structures. Nevertheless, contrary to the claim by S12, 
the large-scale field topology, characterised by the radial field com- 
ponent, was very similar for the maps obtained with the optimal 
and artificially large regularisation, so the Tikhonov regularisation 
operated as expected. KW10 have made it clear that the two values 
of the regularisation parameter are by no means equivalent because 
applying the higher regularisation value yields a far worse fit to the 
observational data. S12 seem to have ignored these explanations 
and the context of KW10 inversions with different regularisation 
parameters, leading to their misinterpretation of the MDI results. 

2.3 Line profile modelling and interpretation of DI maps 

The general goal of Doppler imaging - reconstructing a two- 
dimensional stellar surface map - is always achieved through a 
modelling of the rotational modulation of spectral line profiles. 
However, it is important to discriminate between the tasks of ob- 
taining a map of a certain parameter characterising the stellar sur- 
face and interpreting this map in terms of a physical quantity. In 
the early days of chemical spot mapping of Ap stars, DI stud- 
ies reconstructe d local equivalent width s instead of real elemen- 
tal abundances ( Goncharskii et al. 1983). The majority of current 
DI and ZDI studies of cool active stars map brightness instead 
of temperature, using a coarse temperature-i ndependent analyti- 
cal description of th e local line profiles (e.g., iDonati et "all l2003l : 
iMarsd en et al. 2 0lTI). Modern abundance DI studies of Ap stars 
fe.g jKochukhov et al.l2004bl : iRice. Holmgren & Bohlende3l2004l : 
iLiiftinger et al.l2010ah are normally based on the theoretical spectra 
calculated for different chemical composition using the same mean 
model atmosphere. The resulting maps are presented in the form 
of horizontal distributions of chemical abundance. But it is im- 
portant to remember that even in this case "elemental abundance" 
is essentially a parameter controlling the local line strength and, 
therefore, it is the former stellar surface property that is mapped 
by an inversion code. Interpretation of the line strength parameter 
can be viewed as a separate analysis step, largely decoupled from 
the mapping itself. Assumptions about the atmospheric structure of 
the regions characterised by different metallicity made in this lat- 
ter step are not necessarily undermining reconstruction of the two- 
dimensional maps and have no possibility of changing the topology 
of the main stellar surface features as long as the resulting local pro- 
files can match those produced with a physically correct and com- 
plete model by introducing an offset of chemical abundance from 
its true value. 

If a magnetic field is reconstructed simultaneously with chem- 
ical maps, there is a similar decoupling between polarisation profile 



4 O. Kochukhov, G. A. Wade and D. Shulyak 



modelling and interpretation of the Stokes / line strength. Even if 
the latter is in error by a significant amount, any decent MDI code 
would compensate for this problem by arriving to a different lo- 
cal abundance, which yields the local polarisation profiles closely 
matching the true ones for the same magnetic field strength and 
orientation. One the other hand, an unaccounted variation of the 
continuum brightness may lead to a systematic over or underesti- 
mation of the magnetic field strength in bright, respectively dark, 
surface regions. 

This robustness of DI with respect to errors in interpretation of 
the local line profiles and decoupling between magnetic and abun- 
dance mapping was ignored in the discussion by SI 2. Throughout 
their paper these authors incorrectly imply that any difference be- 
tween the local profiles corresponding to an approximate versus 
a fully realistic treatment of the line formation will automatically 
translate to major surface structure artefacts, invalidating magnetic 
mapping. This viewpoint is fundamentally flawed because it ne- 
glects to take into account the intrinsic ability of an MDI code 
to compensate the line intensity errors arising from the use of an 
approximate relation between the local line strength and chemical 
composition. 

2.4 An evidence for small-scale fields in Ap stars 

Traditionally, analyses of magnetic field topologies of early-type 
stars relied on fitting low-order multipolar field models to the phase 
curves of integral observables, such as the mean longitudinal field 
and the mean field modulus, inferred from th e low- or moderate- 
resolution Stokes I and V observations dLandstreet & Matrivsl 
2000; Bag nulo et alJl2002h . Extension of these studies to broad- 
band linear polarisation revealed significant deficiencies of the mul- 
tipolar modelling approach, hinting at the presence o f more com- 
plex magnetic topologies on the surfaces of Ap stars (iLerov et alj 
Il995l : ILerov. Landolfi & Landi degl 'Innocent! Il996h . A direct 
comparison of the multipolar field geometry predictions with the 
high-resol ution MuSiCoS four Stokes parameter observations of 
Ap stars (IWade et al.ll2000l : iBagnulo et alJbOOlh demonstrated a 
major failure of these topologically simple models in reproducing 
the shapes and amplitudes of the observed Stokes Q and U profiles. 

Inversions based on these pioneering full Stokes vector ob- 
servations showed that small-scale field structures superimposed 
on the global bipolar background are re quired to match the linear 
polarisation inside spectral line profiles ( Kochuk hov et al.l l2004a; 
KW10). The KW10 analysis of a 2 CVn provided the most clear-cut 
example of this situation. The MDI based on the Stokes I and V 
subset of the four Stokes parameter observations yielded a dipolar- 
like field geometry. On the other hand, extending the inversion to 
all four Stokes parameters revealed a pair of high-contrast mag- 
netic spots superimposed on the dipolar-like background. As dis- 
cussed above, this is definitely not a regularisation artefact as the 
difference between the linear polarisation profiles corresponding to 
the magnetic models with and without the small-scale features is 
highly significant. 

The claim by S 12 that "a mere 3 out of 20 phases are responsi- 
ble" for this complexity of a 2 CVn magnetic maps is false. KW10 
have unambiguously stated in the text of their paper and clearly il- 
lustrated with Fig. 5 that the discrepancy between the "simple" and 
"complex" magnetic field geometries is large in comparison to the 
observational noise over at least half of the stellar rotational cy- 
cle, corresponding to approximately 10 phases. The model lacking 
small magnetic spots yields a substantially higher chi- square and 
cannot even qualitatively reproduce the double-wave shape of the 



Q and U profiles observed for the 0.3-0.6 phase interval. Below we 
show that for Stokes QUV this difference is actually considerably 
more important than the impact of using individual local model at- 
mospheres, as urged by S12. 

S12 objected the very possibility of resolving small-scale 
surface features in relatively slowly-rotating stars like a 2 CVn 
(v e sin i = 18.5 kms -1 ) on the grounds that the difference in 
Doppler shifts between the leading and trailing edges of such 
spots is small in comparison to the instrumental resolution. Ac- 
cording to them, the spatial resolution of the stellar DI images 
is entirely determined by the relation of v e sin i and instrumen- 
tal profile width. Any reader familiar with stellar surface map- 
ping will immediately recognise an obvious flow in this reason- 
ing. S12 consider the surface resolution provided by a single ob- 
servation, while DI makes use of many time-resolved spectra (20 
in the case of KW10 study of a 2 CVn), extracting information 
not only from the Doppler shifts but also from rotational modu- 
lation. The latter is more important for magnetic mapping since 
polarisation profiles can exhibit a huge rotational variability even if 
v e sin i is negligible. Numerical experiments by KP02 have demon- 
strated that this allows an application of MDI to all types of early- 
type magnetic stars, regardless of their projected rotational ve- 
locities. Indeed, magnetic DI maps have been obtained for many 
very sl o wly rotating (v e sin i < 5 km s ~ 1 ) early- type ( Donati et all 
2006b; Lufting er et all l2010bh an d late-type dPetit et all 120081 : 
Morin et al. 2008; Fares et al. 2010) stars. A reader may also note 
that with systematic application of the reasoning of SI 2, one is 
bound to conclude that stellar surface mapping is completely un- 
feasible with photometric data - a ludicrous assertion given the 
well- known success of many photometric star spot investigations 
(e.gjKorhonen. Ber dvugina & Tuominen 20021; lLanza et alj 2009: 
iMosser et aDl2009l ; iLiiftinger et al.ll2010ah . For insta nce, the high 
photom etric stability of the CoRoT satellite allowed IMosser et al] 
(2009) to detect photometric signatures of spots as small as a few 
degrees across, whereas according to S12 such a study should not 
have provided any surface spatial resolution at all. 

Summarising, the general arguments by S12 against the de- 
tectability of the complex magnetic fields on Ap stars are demon- 
strably incorrect since they neglect some of the very basic DI prin- 
ciples and ignore a large body of contrary evidence in the literature. 



3 MODEL ATMOSPHERES AND LOCAL PROFILES 

3.1 Atmospheric modelling with the LLmodels code 

Any meaningful analysis of the effects of non- solar chemical com- 
position on the atmospheric structure of Ap stars requires an ad- 
vanced opacity sampling model atmosphere code, which allows to 
carry out accurate calculations for arbitrary elemental abundances. 
In this study we use versio n 8.8 of the well-est ablished model atmo- 
sphere code LLmodels dShulvak et al.l2004l) . This program treats 
atomic line opacity by a direct, line-by-line spectrum synthesis and 
is abl e to account for individu al non- solar chemical abundance pat- 
terns dKhan & Shulyakl 2007) and inhomogeneous vertical distri - 
bution of elements kochukho v et alJl2009tt IShulvak et alll2009h . 
The code also includes a provision for treating the effects of a 
strong magnetic field on the line op acity, including Zeeman split- 
ting and polarised radiative tr ansfer jKochukhov. Khan & Shulyakl 
2005: iKhan & Shulvakll2006h . and a contri bution of the Lorent z 
force to the equation of hydrostatic balance dShulvak et al.ll2007h . 
However, as demonstrated by the previous studies, for the moderate 
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Figure 1. Distribution and visibility of the Fe abundance values in the sur- 
face map derived for a 2 CVn by KW10. Upper panel: Fe abundance for the 
visible part of the stellar surface (latitude ^ 60°). The vertical dashed line 
shows the weighted mean abundance, log N-p e /Ntot = —3.79, used for 
the reference model atmosphere. The dashed lines indicate four additional 
Fe abundances for which we also calculated model atmospheres. Lower 
panel: Visibility function (see Section ing illustrating relative contributions 
of the surface zones with different Fe abundances to the disk-integrated 
spectra of a 2 CVn. 

2-4 kG magnetic field of a 2 CVn, Zeeman splitting and polarised 
radiative transfer can be safely neglected in the model atmosphere 
calculations. 

Using the final Fe abundance map derived for a 2 CVn by 
KW10, we determined a mean surface-integrated Fe abundance of 
log Npe/Ntot = —3.79 using a weight function according to the 
Eq. (64) of PK02. The full amplitude of the Fe abundance variation 
over the stellar surface is 4.9 dex, but the largest values in this range 
are represented by a small number of surface elements and make in- 
significant contribution to the disk-integrated Stokes profile spec- 
tra. Figure Q] shows a distribution of the Fe abundance values for 
the visible part of the stellar surface, which corresponds to the lat- 
itude range from —90° to +60° for the inclination angle i — 120° 
adopted for a 2 CVn. Most of the discussion by S 12 is focused on 
model atmospheres computed for extreme Fe overabundances, such 
as log ATpe/iVtot = —2.0, —1.75, and —1.5. However, the reader 
can readily see from Fig. [T] that these high abundances are quite 
rare. In the actual Fe map published by KW10 only 2.5% of the 
visible surface has log Npe/Ntot ^ — 2. 

In fact, the zones with a high Fe concentration are even less 
important for the disk-integrated Stokes spectra because they are 
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Figure 2. Temperature as a function of the continuum optical depth at A = 
5000 A {upper panel) and Rosseland optical depth {lower panel) in the 
model atmospheres of a 2 CVn computed for different Fe abundances. 



mostly located on the stellar hemisphere directed away from the 
observer. The lower panel of Fig. Q] quantifies this statement by 
plotting a visibility function of different Fe abundances, defined as 
the projected areas of the corresponding surface elements summed 
over all 20 rotational phases of the a 2 CVn observations analysed 
by KW10. If we normalise the area under this visibility function 
to 1, the relative contribution of zones with log N^e/Ntot ^ — 2 is 
only 1.1 x 10 -2 . Thus, contrary to the claims by S12, regions of ex- 
treme Fe overabundance contribute negligibly to the MDI analysis 
of a 2 CVn. 

Aiming to present a balanced analysis of the effects of metal- 
licity on the model atmosphere and spectra of a 2 CVn, we con- 
sidered the Fe abundance range from log Np e /Ntot — —2.19 to 
log NY e /Ntot — —5.91, encompassing 90% of the visible stellar 
surface. We have covered this range with five model atmosphere 
calculations, differing only by the adopted Fe abundance. For 
simpli city, the solar abundances from lAsplund. Grevesse & Sauvall 
(2005) are assumed for all other chemical elements. For all model 
atmosphere calculations reported in this paper we used a 72-layer 
depth grid with an equidistant spacing in the log T5000 range be- 
tween —7 and +2. The frequency-dependent quantities were cal- 
culated between 50 A and 100000 A with a constant wavelength 
step of 0.1 A. The resulting sampling at ~ 10 6 frequency points 
yields more precise line opacities than the Atlas 12 calculations at 
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120000 frequencies employed in the study by SI 2. LLmodels it- 
erations were continued until the temperature change was less than 
1 K at all depths. This condition also satisfied the criteria of the flux 
constancy and radiative equilibrium to 0.1-0.3%. S 12 admitted that 
their model calculations were plagued by convergence problems 
for high Fe abundances. We encountered no such problems with 
LLmodels. 

Our computation of the line absorption coefficient was based 
on the ato mic data for ^ 6.6 x 10 7 lines extracted from the VALD 
data base teupkaetalJll999h . In addition to the standard lists of 
transitions between measured energy levels, this includes a new 
m assive predicted line list calcu lated by R. KurucJ]] and discussed 
bv lGrupp. Kurucz & Tanl d20Q9h . These spectral line data are essen- 
tial for accurate opacity calculations in the context of Ap star atmo- 
spheric modelling, especially for large Fe overabundances explored 
in this study. It is not clear what predicted line list was employed 
by SI 2. 

The distribution of temperature with optical depth for five 
LLmodels atmospheres computed for different Fe abundance is 
illustrated in Fig. [2 Looking at the plot of temperature as a function 
of T5000 , we see a prominent backwarming, which leads to a steeper 
temperature gradient and changes the temperature at log T5000 = 
by 500-1500 K. These effects are less extreme in the higher atmo- 
spheric layers corresponding to the typical line-forming regions. 

It is important to realize, however, that T5000 reflects changes 
in the model structure only at one wavelength. Thus, examining the 
plot of T vs. T5000, one can only say that for the models with dif- 
ferent iron abundance the formation depth of the continuum flux at 
A = 5000 A corresponds to substantially different local tempera- 
ture. But this picture and the resulting conclusions may change if 
one chooses a monochromatic optical depth which is weakly af- 
fected by the Fe abundance. 

A more natural way to compare model atmosphere structures 
is to plot physical quantities against an optical depth computed with 
some mean opacity coefficient. In particular, an optical depth based 
on the Rosseland opacity coefficient k-r has a certain advantage: 
spectral regions with the smallest opacities transfer most of the ra- 
diative energy and at the same time contribute most to kr. Exam- 
ining the physical quantities of different models as a function of 
tr allows one to trace changes of the model structure on a natu- 
ral depth scale representing frequency-integrated properties of the 
radiation field. 

The lower panel in Fig. [2] illustrates the depth dependence of 
temperature for the same set of LLMODELS atmospheres as men- 
tioned above, but now as a function of tr. It is obvious that the ef- 
fects of different iron abundance introduce much smaller changes 
compared to those implied by the T vs. T5000 plot. In particular, the 
temperature variation between different models is reduced to only 
a few hundred K. In the sub-photospheric layers T(tr) is very sim- 
ilar for all models simply because of the fundamental property of 
tr: when tr ^> 1 the temperature distribution asymptotically ap- 
proaches that of a grey model, which, as is well known, depends 
only on T e ff and tr. Therefore, regions of the same tr exhibit the 
same temperature. 

Concluding, a huge variation in the model temperature dis- 
tributions plotted against T5000 should not be mistaken for an ex- 
haustive and complete comparison of different models. A complete 
analysis should always include not only comparisons of the depth- 
dependent thermodynamic properties, but also different spectro- 
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scopic and photometric observables (such as energy distributions, 
hydrogen lines, photometric parameters, metallic line spectra) in 
order to quantify the importance of detected changes in the context 
of particular investigations. 

Let us note in passing that modification of the atmospheric 
structure due to enhanced metal abundances is closely connected 
to a well-known effect of the redi stribution of st ellar flux from 
UV to longer wavelength regions jLeclmmd [l973k making the 
Fe spot region darker in the UV and brighter in the optical 
w avelengths than the su rrounding atmosphere (e.g., see Fig. 9 
of ILiiftinger et Contrary to the impression conveyed 

by SI 2, this prominent influence of chemical composition on the 
structure of stellar surface layers is well-known to the Ap star 
community and has been addressed in numerous applications of 



2009b 



2010b; 



the LLmodels code to individual s t ars (j Kochukhov et al] 
IShulvak. Kochukhov & Khanl 120081 : IShulvak et al.1 I2009L 
ILiiftinger et alJuOlOal) - a significant and fertile research direction, 
seemingly unnoticed by SI 2. 



3.2 Local line profiles and continuum intensities 

To investigate repercussions of metallicity-induced changes of at- 
mospheric structure on the local line shapes and continuum in- 
tensity, we computed the Stokes IQUV profiles of the Fell 
5018.4 4 A line using polarised spectrum synthes is code Syn- 
MAST feochukhov. Makaganiuk & Piskunovll2010l) . This spectral 
line choice is determined by previous magnetic DI studies (K04, 
KW10), where this feature, along with Fell 4923.93 A, was one 
of the very few individual metal lines showing detectable lin- 
ear polarisation signature in the MuSiCoS four Stokes param- 
eter spectra of Ap stars JWade et all l200(t) . As emphasised by 
KW10, Fell 4923 and 5018 A generally represent a poor choice 
for Doppler imaging of chemical spots because these lines are 
saturated in a typical spectrum of a Fe-rich Ap star and thus 
are weakly sensitive to the horizontal abundance variations. On 
the other hand, strong li nes are noticeably affected by vertical 
abundance stratification iRvabchiko va. Leone & Kochukhov| [2005; 
Kochukhov et al ] l2006h . which is usually ign ored in D oppler 
mapping. Non-magnetic abund ance DI analyses dKochukhov et alj 
l2004bl : ILiiftinger et al.1 1201 Oah as well as MDI studies based on 
the Sto kes I and V spectra dKochukhov et al.l2002l ; lLuftinger et alj 
2010b) normally avoid such strong lines, dealing with weak or in- 
termediate strength features. To assess how the latter lines behave 
in the Fe-rich LLmodels atmospheres, we have computed a sec- 
ond set of Stokes profiles for a fictitious weaker iron line with the 
same parameters as Fell 5018.44 A, but with log gf reduced by 
1 dex. 

Starting from the non-magnetic case, we show in Fig. [3^ nor- 
malised disk-centre Stokes / profiles of Fell 5018.44 A and the 
weaker line for all five model atmospheres discussed above. One 
set of line profiles was computed self-consistently with the Fe 
abundance adopted for the model atmosphere calculation, while 
another set represents the usual approach based on a mean atmo- 
spheric model with log NF e /N to t = —3.79 and different Fe abun- 
dances in the spectrum synthesis. These calculations were repeated 
for a 3.5 kG longitudinal field as well as for a 4.0 kG transverse 
field in order to assess any impact of the atmospheric structure on 
the Stokes V and Q profiles, respectively. The V and Q spectra 
are illustrated in Fig. [3]d and c. All calculations included convolu- 
tion with a Gaussian instrumental profiles appropriate for the four 
Stokes parameter MuSiCoS observations employed by K04 and 
KW10. 
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Figure 3. Comparison of the normalised disk-centre profiles of the Fe II 5018.44 A line (lower profiles in each panel) computed self-consistently for the model 
atmospheres with different iron abundances {solid lines) with the profiles calculated for the same abundances and a single mean atmosphere {dashed lines). 
The same comparison is presented for a fictitious weaker Fe II line (upper profiles in each panel). The corresponding spectra are shifted vertically by 0.2 of the 
continuum. The three rows of plots show calculations for a) Stokes / in the absence of magnetic field, b) Stokes V for the 3.5 kG line of sight magnetic field, 
c) Stokes Q for the 4.0 kG transverse magnetic field. 
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Figure 4. Maximum relative profile difference of the normalised disk-centre 
Stokes / spectra of the Fe n 5018.44 A line calculated self-consistently with 
a model atmosphere structure and using a single mean model atmosphere. 
Symbols show calculations for zero magnetic field {circles), 3.5 kG line of 
sight field {triangles), and 4.0 kG transverse field {squares). 

Figures [3^-c leave no doubt that the difference between the 
normalised local line profiles computed with different model atmo- 
spheres and with a mean atmosphere and different Fe abundances 
is extremely small. One cannot visually distinguish the two sets of 
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Figure 5. Same as Fig.|4]for the maximum relative profile difference of the 
normalised disk-centre Stokes V in the 3.5 kG line of sight field {circles) 
and Stokes Q in the 4.0 kG transverse field {triangles). 

calculations except for the strongest line, forming in the most Fe- 
rich atmosphere. Even in that case, the discrepancy is only 4% of 
the continuum for the line which is 83% deep. In general, the max- 
imum relative difference does not exceed several % of the Stokes / 
profile amplitude for Fell 5018.44 A (Fig.© and is entirely negli- 
gible for a weaker line. 
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Figure 6. Relative continuum intensity at A = 5018.44 A as a function of 
the Fe abundance adopted for the model atmosphere calculation. 



How does this difference compare with the sensitivity of 
Stokes I to abundance variation? Considering the Fell 5018.44 A 
line for the model with log N^e/Ntot = —2.19, we found that the 
Fe abundance needs to be changed by just « 0.10 dex, or 2% of 
the logarithmic Fe abundance range for a 2 CVn, to obtain the same 
maximum profile difference as seen in Fig. [4] 

The metallicity-induced modification of the atmospheric 
structure leads to a tiny variation of the normalised Stokes V and Q 
profiles, illustrated in Fig. [3]3 and c, respectively. In the worst case 
scenario, these changes amount to about 0.05-0.10 of the full po- 
larisation profile amplitude (see Fig. [5}. Generally, the polarisation 
profiles become slightly broader, while their low-order moments 
remain unchanged. 

The full Stokes profile calculations presented in this section 
fully support the conclusion by KW10 that "the local line profiles 
are sensitive to model structure effects to a much smaller degree 
than to changes of abundance or magnetic field". At the same time, 
enhanced metallicity leaves a major imprint on the continuum in- 
tensity, which is increased in the optical due to the flux redistri- 
bution from the short- wavelength regions. Figure [6] demonstrates 
that the Fe-rich zones are expected to be up to 40% brighter in 
the continuum at A 5018.44 A with respect to the reference model, 
while the regions with a relative Fe underabundance turn out to be 
up to 10% fainter. The consequences of these continuum intensity 
changes for the disk-integrated Stokes IQUV profiles and for the 
reconstruction of the surface abundance and magnetic field distri- 
butions will be addressed in the next section. 

The general conclusion regarding the insignificance of the dif- 
ferential line-blanketing for the local spectral line shapes emerg- 
ing from our calculations seems to be drastically different from the 
claims made by S 12 in Section 4 of their paper. This discrepancy 
comes from the fact that S12 have neglected to consider the en- 
tire Fe abundance range and different line strengths, limiting their 
discussion to the profile of a single, very strong line computed for 
extreme Fe overabundances which occupy a tiny fraction of the 
visible surface of a 2 CVn. In addition, their local profile analysis 
lacked a comparison of the intensity differences with the actual pro- 
file amplitude, as in our Fig. [3] and did not at all address polarisa- 
tion spectra. As revealed by our more comprehensive calculations, 
the atmospheric structure variations result in a systematic change 
of the continuum fluxes but has a minuscule influence on the nor- 
malised intensity and polarisation lines shapes, except for a small 
region of the parameter space (the strongest line in the most Fe-rich 



atmosphere). As discussed above, we do not expect this region to 
contribute significantly to the disk-integrated line profiles. This is 
directly confirmed by the spectrum synthesis calculations presented 
in the next section. 



4 MAGNETIC DI WITH SELF-CONSISTENT MODEL 
ATMOSPHERE STRUCTURE 

4.1 Magnetic inversion code Inversl3 

We have carried out a series of forward line profile calculations to 
assess the impact of the anomalous local model atmosphere struc- 
ture on the disk-integrated Stokes spectra of a 2 CVn. For these ex- 
periments we used our new magnetic DI code INVERS13. This ef- 
ficient and versatile code represents a further development of our 
magnetic inversion tool InversIO , which has been describ ed in 
detail by PK02, evaluated by iKochukhov & Piskunovl d2002h . and 
was subsequently applied in several studies of a non-uniform distri- 
bution of magnetic fields and chemical abundances on the surfaces 
of Ap stars (K04, KWlO. lLuftinger et alJl2010bl : IKochukhov et alJ 
120111) . 

Invers 13 and InversIO adopt similar regularisation and 
non-linear optimisation strategies, based on the Tikhonov regular- 
isation and modified Levenberg-Marquardt method, respectively. 
However, instead of using a single model atmosphere structure as 
was done in all previous DI studies of Ap stars, INVERS 1 3 is capa- 
ble of mapping the three components of the magnetic field vector 
and treating simultaneously one additional scalar surface map (tem- 
perature, abundance of a given element, magnetic field strength) 
fully self-consistently with the local model atmosphere structure. 
The basic methodology of applying INVERS13 to Ap stars is thus 
very similar to a self-consistent magnetic and temperature mapping 
of cool active stars, except that the parameter distinguishing differ- 
ent local model atmospheres is a chemical abundance instead of 

At the start of calculations, a grid of model atmo- 
spheres (either ATLAS, MARCS or LLMODELS) is pro- 
cessed by INVERS 13. The equations of chemical equilib- 
rium (cool atmospheres) or ionisation balance (hot atmo- 
spheres) are solve d with the help of a robust mole cular 
equilibrium solver JValenti. Piskunov & Johns-Krulll [l998) im- 
plemented in SME dValenti & Piskunovl Il996h and in the 
Synth3/Synthmag/Synma st family of codes (iKochukhovl 
120071 : IKochukhov et all l2010h . Then, INVERS13 calculates the 
depth-dependent line and continuum opacity coefficients for all 
model atmospheres and all molecular and atomic spectral lines 
included in th e input list. The a tomic line parameters are taken 
from VALD (Kupka etail Il999h . whereas the data on molec- 
ular transitions are ta ken from the MARCS opacity database 
dGustafsson et al.ll2008h . 

The core module of INVERS 1 3 is the quadratic DELO algo- 
rithm for highly accurate solution of the LTE polarised radiative 
transfer equation in an arbitrary stellar atmosphere permeated by 
a depth-independent magnetic field with a given (local) strength 
and orientation. PK02 demonstrated that the quadratic DELO for- 
mal solver is superior in terms of the speed and accuracy to the 
Zeeman Feautrier method used by SI 2. For each element in a stel- 
lar surface grid optimally sampling different latitudes (see Fig. 5 
in PK02), our code evaluates the line absorption matrix by adding 
opacities of all contributing Zeeman components and computes the 
local Stokes IQUV profiles and continuum intensity I c by per- 
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forming a quadratic interpolation within three sets of model spec- 
tra corresponding to the model atmospheres bracketing the value 
of a scalar parameter. This also yields precise derivatives of the 
Stokes IQUV spectra and associated intensity I c with the respect 
to a scalar parameter distinguishing different models in the input 
model atmosphere grid. One-sided derivatives with respect to the 
radial, meridional, and azimuthal magnetic field vector components 
are computed numerically. The full calculation of the local Stokes 
spectra and their derivatives thus consists of 12 separate polarised 
line profile calculations for every surface element at every relevant 
rotational phase. 

The local profiles are convolved with a Gaussian function 
to model instrumental and radial-tangential macroturbulent broad- 
ening. Profiles are Doppler- shifted and summed with appropriate 
weights for each rotational phase. A similar phase-dependent disk- 
integration is carried out for the continuum intensity. The resulting 
Stokes IQUV spectra are normalised by the respective continuum 
fluxes. The matrix of partial derivatives of every observed spectral 
pixel with respect to every surface parameter is established taking 
into account both the magnetic and model structure sensitivity of 
the line profiles and the model atmosphere sensitivity of the contin- 
uum fluxes. Using these derivatives, the code establishes the curva- 
ture matrix required by the Levenberg-Marquardt algorithm, adds 
regularisation terms and iteratively adjusts surface maps to match 
the observations. 

INVERS13 takes advantage of the standard MPI libraries 
available at any modern multi-CPU computer to perform highly 
efficient parallel radiative transfer calculations with a novel master- 
slave technique discussed by PK02. We usually run the code with 
one master process collecting results of the radiative transfer calcu- 
lations by 16-32 slave processes. The code was successfully tested, 
showing a linear scaling, with up to several hundred CPUs in simul- 
taneous inversions of over a thousand atomic and molecular lines 
spread over ten wavelength intervals. In practice, the number of 
spectral lines and wavelength intervals modelled simultaneously is 
limited only by the available RAM. For the calculations described 
here, modelling of one or two Fe II spectral lines at 20 rotational 
phases required less than one minute of computing time for the 
forward spectrum synthesis and about three hours for a complete 
inversion on a modest, 10-year-old 16-CPU computer. This experi- 
ence shows that currently available computing resources and stan- 
dard scientific programming techniques based on well-established 
Fortran codes enable a fully realistic treatment of the lateral vari- 
ation of stellar atmosphere, without a recourse to the next gener- 
ation computers or exotic programming languages advocated by 
S12. 

INVERS 1 3 was rigourously tested with numerical experiments 
for cool active star s using simulated circular pola risation and full 
Stokes vector data dKochukhov & Piskunovl 120090 . The code was 
also applied to interpret circular polarisation obs ervations of the 
active RS CVn star II Peg ( Kochuk hov et al.l20 09a). A comprehen- 
sive MDI study of this target is currently underway. 

4.2 Disk-integrated IQUV profiles of the Fe n 5018 A line 

We have used the forward polarised spectrum synthesis calcula- 
tions with INVERS10 and INVERS13 to directly answer the ques- 
tion of how important is the local atmospheric structure for the 
KW10 analysis of the Fe abundance distribution and magnetic 
field geometry of a 2 CVn. We adopted the magnetic and abun- 
dance maps derived in that study and employed a single Atlas 9 
model atmosphere with parameters T e ff = 1 1600 K, log g = 3.9 and 
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Figure 7. Theoretical four Stokes parameter profiles of the Fe II 5018.44 A 
line corresponding to the final magnetic and abundance maps derived for 
a 2 CVn by KW10. This plot compares the self-consistent profiles (solid 
line) with the standard calculations using a mean model atmosphere (dashed 
line). For comparison, we also show profiles for the smoothed magnetic 
field geometry model discussed by KW10 (dotted line). Spectra for consec- 
utive rotational phases are shifted vertically. Rotational phases are indicated 
on the right. The bars at the lower right of each column show the horizon- 
tal and vertical scale (0.5 A and 1% of the Stokes I continuum intensity 
respectively). 
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Figure 8. Phase dependence of the amplitude, observational uncertainty, 
and the difference between theoretical Stokes parameter profiles illustrated 
in Fig. [7] This plot compares the amplitude of Stokes profiles (black solid 
line) with the discrepancy between the standard and self-consistent calcu- 
lation (dashed line), and observational uncertainty (dotted line). For com- 
parison, we also show the difference between the standard profiles and the 
spectra corresponding to the smoothed model of the magnetic field topology 
(red solid line). 



[M/H] — +1.0 for the full Stokes parameter synthesis of the Fe II 
5018.44 A line with InversIO. These calculations are identical to 
the final Stokes profile fits presented by KW10. In particular, we 
consider the same set of rotational phases and take into account the 
spectral resolution of the MuSiCoS Stokes vector data available for 
a 2 CVn. 

Another set of line profiles was produced with the Invers 1 3 
code using the grid of LLmodels atmospheres discussed in Sec- 
tion [3j] To avoid extrapolation in Fe abundance, we augmented 
this grid with two extreme models computed for log N^e/Ntot — 
— 1.60 and —6.40. We checked that using a denser grid of 13 
model atmospheres makes no difference. The resulting normalised 
Invers 1 and Invers 1 3 Stokes parameter profiles are compared 
in Fig. |7] In addition, we plot a set of synthetic InversIO Stokes 
IQUV parameters corresponding to the smoothed model of the 
magnetic field geometry of a 2 CVn discussed by KW10. A de- 
tailed assessment of the difference between the "standard", "self- 
consistent" and "smoothed" model profiles with respect to the ob- 
servational uncertainties is presented in Fig. [8] We note that for 
polarisation this comparison represents a conservative estimate of 
possible errors because we do not allow the code to adjust the 
Stokes / profile intensity as it would do in a real magnetic inversion 
(see discussion in Section [23]), thereby improving the fit to QUV. 

What can we learn from this comparison? Clearly, the effect of 
self-consistent calculations is noticeable for Stokes I. The intensity 
profile changes by 2-4% of the continuum, which is not negligi- 
ble compared to the central line intensity of 27-35%. In agreement 
with the prediction made in Section [3T2l the self-consistent inten- 
sity profiles are systematically deeper due to an increased contin- 
uum brightness of the Fe-rich surface regions. However, this effect 



can be compensated by a mere ~ 0. 1 dex reduction of the average 
Fe abundance, so it is not immediately obvious that accounting for 
the local atmospheric structure should lead to any major change of 
the abundance inversion results and remove extremely high abun- 
dances from the Fe map. 

At the same time, it is striking to see how little effect the use 
of self-consistent modelling of the Stokes profiles has on the linear 
and circular polarisation spectra. The discrepancy in Stokes V is 
only 0.2-0.4%, which exceeds the observational uncertainties but 
is negligible compared to the maximum circular polarisation signal 
of 9%. The impact on the linear polarisation spectra is even smaller. 
It is evident from Fig. [8] that accounting for the local atmospheric 
structure changes Stokes Q and U by about 0.1%, which is below 
the noise level of observations for all but two rotational phases. 

It is enlightening to compare the difference between the stan- 
dard InversIO calculation and self-consistent Invers 13 model 
with the discrepancy between the former spectra and those corre- 
sponding to the smoothed version of the magnetic field topology 
of a 2 CVn (thick solid line in Fig. [8]). We remind the reader that 
this discrepancy, constituting the main evidence for the small-scale 
magnetic fields in a 2 CVn, was dismissed by S 12 as insignificant. 
Figures [7] and [8] leave no doubt that the conclusions of S12 are er- 
roneous. A signature of the small-scale fields yields a systematic 
modulation of the linear polarisation spectra, which is above the 
observational noise for 15 out of 20 rotational phases in Stokes 
Q and for 18 phases in U. We see that the actual situation is just 
the opposite of that claimed by S12 in Sections 5 and 7 of their 
work: the influence of structured fields is highly significant for the 
Stokes VQU profiles, while the assumption of a mean model atmo- 
sphere is largely unimportant for V and completely negligible for 
QU even in this pessimistic assessment (no abundance adjustment, 
a strong line). 

The results of the forward full Stokes profile calculations un- 
equivocally demonstrate that the speculations by S12 about the sig- 
nificance of the local model atmosphere structure for the modelling 
of polarisation in a 2 CVn are incorrect. At this point the reader is 
certainly entitled to ask what led S12 to these erroneous conclu- 
sions. The Fe distribution considered by S12 in Section 5 of their 
paper appears to be far more extreme (log(A^Fe/^V to t) ^ — 2 over 
18% of the visible stellar surface right at the disk centre) compared 
to the one found for a 2 CVn by KW10 (only 2.5% of the visible 
surface close to the limb). The unspecified "eccentric dipole" with 
an unusually high contrast (field strength variation by a factor of 
3.7 instead of a factor of 2 expected for a centred dipole) seems to 
have been hand-picked to maximise profile differences. Curiously, 
describing the Stokes parameter calculations for this large Fe spot 
and unusual magnetic geometry, S12 do not mention taking into 
account the instrumental profile of the spectrograph, although they 
stress a moderate resolution of the MuSiCoS four Stokes parame- 
ter data elsewhere in the paper. It is also very surprising that the 
Q and U profiles, for which the effect of self-consistent modelling 
is the least important but which conversely have played a key role 
in the discovery of small-scale magnetic features on the surface of 
a 2 CVn, are completely neglected in that part of the S 12 paper. 

Calculations presented in Section 7 of the S 12 paper are sim- 
ilarly inadequate. In that part of their study the authors claimed 
to have assessed the feasibility of detecting a representative sig- 
nal of a high-contrast magnetic spot. However, the scenario consid- 
ered is not applicable to the case of a 2 CVn because the adopted 
spot diameter is underestimated by a factor of two and the im- 
posed correlation of the small-scale abundance and magnetic field 
enhancements is contrary to what was found for that star. Further- 
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Figure 9. Surface distribution of Fe abundance derived for a 2 CVn from the Stokes IQUV profiles of the Fe II 4923.93 and 5018.44 A lines, a) Inversion 
with InversIO using a single model atmosphere (KW10), b) inversion with Invers13 using a self-consistent treatment of the Fe abundance and model 
atmosphere structure (this study), c) the difference map. The spherical maps are shown for five equidistant rotational phases as indicated on the top of the 
figure. The aspect corresponds to the inclination angle i = 120° and vertically oriented rotational axis. The contour lines are plotted with a step of 0.5 dex. 
The stellar rotational equator is shown with a thick line, while the vertical bar indicates rotational axis. 



more, S12 produced a small-scale magnetic feature by simply scal- 
ing all three vector components in an arbitrarily chosen location 
on the star instead of introducing horizontal field structures similar 
to those found for a 2 CVn by KW10. These pitfalls notwithstand- 
ing, the authors obtained differential polarisation signatures of up to 
^0.2% (for a spectral resolution of R = 35000). These signatures 
are about twice lower than in the real star (see Q and U panels in 
Fig. [8}, yet could still be detected with the help of multiple spectral 
lines and multiple Stokes profile observations at S/N ~ 10 3 . But 
this profoundly multi-line, multi-phase character of modern MDI 
was not considered by SI 2, who limited their analysis to a single 
spectral line at one rotational phase. 

In summary, our forward Stokes profile calculations demon- 
strate a moderate metallicity effect of the local model atmosphere 
on the Stokes / profile of one of the strongest Fe lines in the optical 
spectrum of a 2 CVn, but very small (V) or negligible (QU) influ- 
ence on other Stokes parameters. Local Stokes parameter calcula- 
tions (Section lX2l l suggest that the corresponding effect on weaker 
lines will be negligible for all four Stokes parameters. 



4.3 Self-consistent inversion of the Fe abundance and 
magnetic field structure 

Taking advantage of the capabilities of Invers 1 3, we have carried 
out a complete reconstruction of the Fe spot distribution and mag- 
netic field geometry on the surface of a 2 CVn using the full Stokes 
vector observations by KW10. Incidentally, this is the first ever at- 
tempt to directly incorporate variation of the stellar atmospheric 
structure in the Doppler imaging of an Ap star. We modelled si- 
multaneously the Stokes IQUV profiles of the Fell 4923.93 and 
5018.44 A lines, keeping the same strength of the Tikhonov reg- 



ularisation and adopting the same stellar parameters as employed 
by KW10 for their final full Stokes vector MDI inversion. The 
iron abundance was treated self-consistently with the model at- 
mosphere structure using the LLmodels grid described above. A 
comparison of the resulting DI map with an outcome of the pre- 
vious INVERS 1 modelling based on a single mean model atmo- 
sphere provides the most direct and conclusive assessment of the 
possible adverse influence of the usual simplified treatment of the 
Ap-star atmospheric structure. We note that the claims of S12 are 
not support ed by a similar assessment because their inversion code 
Jstif3l996) is not capable of magnetic mapping and, at the moment, 
cannot account for the lateral variation of atmospheric structure. 

A comparison of the Fe abundance maps reconstructed by 
INVERS 10 and INVERS 13 from the same observational material 
is presented in Fig. [9] Clearly, the chemical spot distributions re- 
covered with the two methods are very similar and there is no sign 
whatsoever of a major re- arrangement of spots or significant re- 
duction of the abundance contrast anticipated by S12. The main 
Fe features are located in the same regions of the stellar surface 
and show consistent amplitudes. We have to turn our attention to a 
difference map (Fig. |9j:) for a closer inspection of the discrepancy 
between the two inversions. Generally, INVERS13 yields a some- 
what lower Fe abundance with respect to INVERS 1 0. The mean dif- 
ference is —0.12 dex, while the standard deviation is 0.33 dex. As 
evident from the difference map, the discrepancy reaches ^0.5 dex 
over a significant fraction of the visible stellar surface. The extreme 
outliers, occupying a few surface pixels, are showing abundance 
differences at the level of ±0.7 dex. Although this change of the 
chemical map is not negligible, it still corresponds to only 10-15% 
of the full logarithmic range of the Fe abundance. At the same time, 
the mean abundance changes by just 2% of that range. 
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Figure 10. Rectangular maps of the radial, meridional, and azimuthal vector components of the magnetic field structure derived for a 2 CVn from the Stokes 
IQUV profiles of the Fen 4923.93 and 5018.44 A lines, a) Magnetic inversion with InversIO using a single model atmosphere (KW10), b) magnetic 
inversion with Invers 1 3 using a self-consistent treatment of the Fe abundance and model atmosphere structure (this study), c) the difference map. The dotted 
line shows the highest visible latitude for the inclination angle i = 120°. The contours of equal field strength are plotted with a 0.5 kG step. 



The application of INVERS 13 to a 2 CVn clearly demon- 
strates that a less detailed treatment of the atmospheric structure by 
InversIO does not introduce spurious abundance structures and 
recovers basically the same Fe horizontal abundance contrast as a 
more sophisticated modelling. We are far from insisting that these 
Fe spot mapping results are completely free from systematic errors. 
But it is obvious that a mean atmosphere assumption is neither a 
dominant nor a significant source of uncertainty for the published 
MDI studies of a 2 CVn and other Ap stars. In particular, this as- 
sumption is definitely not the reason for the extremely large over- 
abundances inferred to exist in the atmospheres of Ap stars, which 
S12 found so difficult to believe. 

The results of the magnetic mapping carried out simultane- 
ously with the reconstruction of the Fe abundance are presented 
in Fig. [T0l It is reassuring to see very little difference between 
the InversIO and Invers 13 magnetic maps, despite a more 
rigourous treatment of the stellar atmosphere by the latter code. 
Even the difference map (Fig. [TUb ) hardly shows any significant 
structures. The mean difference is about 0.05-0.10 kG for all three 
magnetic field vector components, while the standard deviation is 
0.19, 0.24 and 0.13 kG for the radial, meridional and azimuthal 
field, respectively. This corresponds to a typical local magnetic 
field error of only 2-3% of the full field amplitude, while even the 



largest difference of ^0.5 kG yields a 7% local field strength er- 
ror. The fine structure of the stellar magnetic field, initially found 
using Invers 1 and most clearly seen in the rectangular map of 
the azimuthal field component, is also faithfully reproduced by 
Invers 13. 



This excellent agreement between the magnetic maps recov- 
ered using different approaches to the atmospheric structure proves 
that MDI based on high-quality linear and circular polarisation data 
is remarkably robust and is insensitive to the possible errors in in- 
terpretation of the Stokes I spectra related to using a single model 
atmosphere. Our results convincingly demonstrate that the bleak 
S12 assessment of this magnetic inversion procedure is grossly in 
error. It appears that their allegations are based on an unreasonable 
extrapolation of a very limited set of forward Stokes parameter cal- 
culations for an artificial abundance distribution not resembling any 
real Ap star. The outcome of the inversions proper, presented here, 
proves that the methodology used by S 12 is misleading and erro- 
neous. A general conclusion, which inevitably follows from this 
discussion, is that the S12 assertions are largely irrelevant in the 
contex t of modern MDI studies by K WlO jKochukhov et all (l200l 
l2004ah . and lLuftinger et all (l2010bh . 
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5 CONCLUSIONS AND DISCUSSION 
5.1 Summary of main results 

In this paper we have confronted the criticism by Stift et 
S 12) of the basic methodological foundations and specific results of 
magnetic and abundance Doppler imaging (MDI, DI) of Ap stars. 
We demonstrated that the doubts about general issues, such as sta- 
bility of magnetic inversions, robustness of the regularisation pro- 
cedure, and detectability of small-scale stellar surface features are 
unwarranted because they were adequately addressed by a number 
of previous publications ignored by S12. Many successful appli- 
cations of MDI to slowly rotating stars and the very possibility of 
photometric mapping of star spots show that the theoretical esti- 
mate by S 12 of the surface resolution attainable by stellar surface 
imaging is wrong. 

For the specific case of the MDI analysis of a 2 CVn by KW10, 
we showed that S12 have greatly overestimated the contribution of 
the surface zones with a high Fe abundance to the disk-integrated 
Stokes spectra and have thus performed their model calculations 
for an unrealistic parameter range. Using a state-of-the-art opacity 
sampling model atmosphere code and considering representative 
Fe abundances, we demonstrated that the effect of metallicity de- 
pendence of the model atmosphere structure is negligible for all lo- 
cal Stokes profiles except those corresponding to the strongest lines 
formed in the most Fe-rich regions (log(ivVe/iVtot) ^ —2). At the 
same time, the flux redistribution characteristic of the Fe-rich mod- 
els leads to a significant continuum brightening of the zones cov- 
ered by overabundance spots - an effect potentially important for 
the spectrum synthesis of Ap stars and for Doppler inversions. 

Using a new magnetic DI code which fully incorporates ap- 
propriate local model atmospheres in polarised spectrum synthesis, 
we performed forward calculations of the four Stokes parameter 
profiles of the Fell 5018.44 A line in the spectrum of a 2 CVn. 
This experiment showed that using individual model atmospheres 
for overabundance regions compared to the usual single model at- 
mosphere assumption yields a small modification of Stokes / but 
is completely negligible for polarisation profiles. In particular, the 
resulting changes of the Stokes Q and U spectra of a 2 CVn do 
not exceed the observational noise of the data used by KW10 and 
are significantly smaller than the signature of small-scale magnetic 
field spots. 

Taking advantage of the capabilities of the new magnetic DI 
code, we have performed the first simultaneous magnetic and Fe 
abundance inversion using a grid of model atmospheres appropri- 
ate for the regions with different elemental abundance. The result- 
ing revised map of the Fe distribution over the surface of a 2 CVn 
changes slightly with respect to the results published by KW10 but 
maintains the same overall spot topology and still shows the very 
high Fe overabundances objected to by S12. Our new inversions 
convincingly demonstrate that, although possibly an artefact of us- 
ing very strong lines, the most extreme Fe concentrations are defi- 
nitely not a result of adopting a single mean model atmosphere. 

At the same time, the magnetic map produced using the new 
inversion method turned out to be practically identical to the one 
inferred by KW10 using a less sophisticated approach. In partic- 
ular, the presence of small-scale magnetic spots on the surface of 
a 2 CVn is confirmed. The new inversions demonstrate that recon- 
struction of the magnetic field configuration is sufficiently decou- 
pled from the abundance mapping to be insensitive to the details 
of model atmosphere structure. These results strongly suggest that 
self-consistent MDI incorporating model atmospheres calculated 
for individual local elemental abundances is an unnecessary lux- 



ury, not justified by the quality of currently available full Stokes 
parameter data. 

5.2 Comparison with ZDI of cool stars 

It is instructive to compare the impact of using a single model at- 
mosphere in the MDI of Ap stars with the repercussions of usual 
assumptions of cool-star ZDI. Although not mentioned explicitly 
by S12, their criticism obviously applies to the magnetic map- 
ping studi es of cool active stars as well. With a few recent ex- 
ceptions ( Carroll et al. 2007; Kochukhov et al. 2009a), these stud- 
ies interpret circular polarisation completely neglecting tempera- 
ture or brightness inhomogeneities. They also adopt a simplified, 
temperature-independent analytical descriptions of the local Stokes 
/ and V line profiles instead of calculations based on realistic 
model atmospheres for different effective temperatures. In contrast, 
before the present study, MDI analyses of Ap stars have used model 
atmospheres to calculate the local Stokes profiles as a function of 
chemical abundance but have ignored continuum intensity varia- 
tions. 

How important is the continuum brightness effect due to 
chemical inhomogeneities compared to that of cool spots on the 
surfa ces of solar-type stars? Using the Marcs model atmosphere 
grid ( Gustafsson et al . 2008), we determined that a 40% contin- 
uum intensity contrast (inferred for a rather extreme case of the 
Fe spots on a 2 CVn) corresponds to « 600 K temperature differ- 
ence between a spot and photosphere for stellar parameters close 
to solar values. This is a moderate temperature contrast relative to 
^phot _ T ST) ot — 1500-2000 K inferred for early-G active stars 
(Berdyugina 2005). This comparison shows that the errors due 
to neglecting continuum intensity variation in the magnetic DI of 
early-type stars are significantly smaller than in the ZDI of cool 
active stars. 



5.3 DI maps and photometric variability of Ap stars 

DI/MDI techniques have been successfully applied to explain ro- 
tational variability of line profiles in various CP stars. The re- 
sulting abundance maps allowed to explore the distribution of 
different chemical elements over stellar surfaces and compare 
these distributions with the predictions of atomic diffusion the- 
ory jMichaud. Charland & Meg essier 1981). On the other hand, for 
many decades magnetic CP stars were known to exhibit photomet- 
ric variations with the same period as the magnetic and spectro- 
scopic changes, but this variability had been never explained quan- 
titatively. The fact that the photometric variability occurs in anti- 
phase in the far-UV and optical spectral regions jMolnad[l97l 
Il975l) suggested that the global flux redistribution caused by the 
phase-dependent line absorption related to chemical spots is the 
most likely cause of the observed photometric variability. But it 
was only recently t hat this hypothesis was rigourously tested by 
iKrticka et al.l (|2007), who presented a successful attempt to inter- 
pret uvby light curves of the hot Bp star HP 3 7776 using DI maps 
of He and Si derived bv lKhokhlova et al.l {2000). The adopted pho- 
tometric variability analysis method was straightforward and real- 
istic: they computed a grid of model atmospheres taking into ac- 
count the local abundances provided by DI maps and then per- 
formed surface integration of the local continuum intensities to 
obtain phase-dependent fluxes in wavelength regions covered by 
different p hotometric filters. T his local model atmosphere analy- 
sis allowed Krticka et al. (120071) to reproduce the amplitudes of the 
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observed light variability in different filters and even to match the 
light curve shapes (which evidently depend on the geometrical dis- 
tribution of the bright and dark areas on stellar surface). 

This work has provided an important and completely indepen - 
dent confirmation of the results of DI bv lKhokhlova etalJ |2000). 
Subsequently, the same procedure of explaining the observed light 
variability by an inhomogeneous surface distribution of chemical 
elements was s uccessfully applied to severa l other magnetic C P 
stars: HR72 24 (IKrticka et al .120091) . g UMa dShulvak et alfcoiOah . 
and CU Vir (IKrticka et al.ll2012h . In all these investigations, the In- 
VERS family of Doppler imaging codes was used to derive abun- 
dance distributions assuming a single model atmosphere in map- 
ping. Furthermore, using the s ame chemical spot imaging method- 
ology, iLuftinger et alj j2010ah have demonstrated a perfect coinci- 
dence between the location of spots of iron-group elements on the 
surface of the Ap star HD 50773 and the distribution of bright re- 
gions inferred from the high-precision C0R0T light curve of this 
star. 

These results obtained from light curve modelling studies def- 
initely proved that an inhomogeneous distribution of chemical ele- 
ments is the dominant source of the observed photometric variabil- 
ity and that the usual DI inversions correctly predict the location 
and amplitude of the main abundance features which give rise to 
photometric variability. Obviously, a major geometrical artefact in 
the spot distribution or a gross overestimate of the Fe-peak abun- 
dances, as suspected by SI 2, would not allow reproduction of the 
stellar light curves. It is quite surprising that this entire series of 
publications, providing a convincing independent validation of the 
abundance DI results, has escaped the attention of SI 2. 

5.4 Final comments 

Besides addressing specific scientific concerns voiced by S12 using 
the new Stokes profile calculations and drawing readers' attention 
to copious scientific literature supporting our conclusions, we are 
led to comment on a number of cases where the S 12 interpretation 
of published studies is subjective and could potentially be mislead- 
ing. 

The choice by S 12 of the abundance DI studies yielding "un- 
realistic" elemental overabundances is highly selective, confined 
to two extreme result s , both published in c onference proceedings 
dKuschnig etal.| [l998: Pisk unov et allll998h . This ignores the ma- 
jority of DI studies published in the refereed literature which re- 
ported much more modest abundance contrasts. 

S12 compare the MDI studies of a 2 CVn bv lKochukhov et alj 
(2002) and KW10 at great length. Excessive emphasis is given to 
the most extreme values of the magnetic and abundance maps; how- 
ever, the maps' information content extends far beyond these ex- 
trema. In fact, the surface distributions agree reasonably well and 
any differences are well understood. KW 10 presented two magnetic 
maps: one obtained from the Stokes IV spectra and another derived 
from all four Stokes parameters. There is no major difference be- 
tw een the former and the r esults of circular polarisation modelling 
bv lKochukhov et all J2002I) . except for a complementary choice of 
the inclination angle by KW10 (Stokes IV data do not allow dis- 
tinguishing i and 180° — i, see, e.g., Bagn ulo et all J2OOO)) and 
a systematically weaker field obtained in the more recent study. 
The latter effect is probably d ue to an excessive longit udinal field 
curve regularisation applied by Kochukh ov et al.l d2002h . The over- 
all global field topology, characterised by the distribution of the 
radial field component, is very similar in all magnetic maps recon- 
structed for a 2 CVn. On the other hand, the discrepancy between 



the two and four Stokes parameter inversions is the central point 
of the study by KW10. This difference indicates a rich information 
content of linear polarisation, not an inve rsion artefact. 

The Fe and Cr maps presented by Kochukhov et al ] d2002h 
and KW10 agree qualitatively. All maps show a relative under- 
abundance around phase 0, with spo ts located on both sides of the 
stellar disk yielding profile doubling. Kochuk hov et al.l J2002) used 
higher resolution data and weaker lines, so their maps are more re- 
liable. The limitations of the abundance mapping using very strong 
metal lines (weakly sensitive to abundance variation and likely af- 
fected by vertical stratification) have been explicitly acknowledged 
by KW10 (Section 4.2): "we do not expect to achieve the same 
quality of abundance inversions as demonstrated by Kochukhov et 
al (2002)". 

Finally, S12 state that the large overabundances inferred from 
DI are at odds with the predictions of theoretical diffusion calcula- 
tions, which have not been shown to produce such large accumu- 
lations. In fact, both cited theoretical diffusion studies have artifi- 
cially limited the maxim um accumulation of chemical elements to 
log(A WA^tot) -2.3 EeBlanc et al.l l2009h or \og(N el /N to t) = 
—3.0 jAlecian & StifJ 120071) in order to avoid numerical insta- 
bilities. However, in both studies this limit has been reached for 
certain chemical ele ments at some optical depths (see F ig. 3 of 
lAlecian & Stif3 d2007h and Fig. 4 of iLeBlanc etaD (l2009h ). Thus, 
the lack of very high elemental overabundances in the equilibrium 
diffusion calculations is an artefact of the limited computational 
capabilities of the currently used atomic diffusion codes. 
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